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The neutron-star inner crust is assumed to be superfluid at relevant temperatures. The 
contribution of neutron quasiparticles to thermodynamic and transport properties of the 
crust is therefore strongly suppressed by the pairing gap. Nevertheless, the neutron gas 
still has low-energy excitations, namely long-wavelength collective modes. We summarize 
different approaches to describe the collective modes in the crystalline phases of the inner 
crust and present an improved model for the description of the collective modes in the 
pasta phases within superfluid hydrodynamics. 
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1. Introduction 

A neutron star has a very dense core, consisting probably of homogeneous and 
very neutron-rich matter (other constituents are protons, electrons, and perhaps 
hyperons; the extremely dense center of the core might even contain deconfined 
quark matter), which is surrounded by the inner and the outer crustPln the crust, 
the baryon density ub is below ~ 0.08 fm” 3 (i.e., mass density g < 1.3x 10 14 g/cnr 3 ) 
and the matter is inhomogeneous, containing dense positively charged “clusters” 
and a highly degenerate electron gas. The difference between the inner and the 
outer crust is that in the outer crust, the clusters are simply neutron-rich nuclei, 
while in the inner crust (at mass density g > 4 x 10 11 g/cm 3 P i.e., baryon density 
ub > 2.4 x 10 -4 fm -3 ), the neutron excess becomes so strong that some neutrons 
are not bound any more in the clusters and form a dilute neutron gas between 
the clusters^ To minimize Coulomb energy, the clusters are believed to arrange in 
a periodic lattice. With increasing density of the neutron gas near the core, the 
crystalline lattice might transform into the so-called “pasta” phases, i.e., clusters 
merge first into rods (“spaghetti”), then into plates (“lasagne”), and then some 
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theories predict also “inverted” geometries where the more dilute neutron gas is 
concentrated in tubes or holes (“Swiss cheese”) before the homogeneous core is 
reached P® 

In the first 50-100 years after the creation of the neutron star in a supernova 
explosion, its core cools down very efficiently by neutrino emission. During this so- 
called crust-thermalization epoch, the crust stays hotter than the core. Since the 
observed temperature is that at the surface, the thermal properties of the crust 
influence the observed cooling curvePln accreting neutron stars, the matter falling 
on the star results from time to time in nuclear reactions in the crust, leading 
to a strong heating (released as X-ray burst). The subsequent cooling of these X- 
ray transients offers another possibility to get information on the thermodynamic 
properties of the crust P 7 

The main ingredients to describe heat transport in the crust are the specific 
heat and the heat conductivity. They depend mainly on excitations whose energy 
is of the same order of magnitude as the temperature T. Since on a nuclear energy 
scale, the temperatures of interest (10 — 100 keV, corresponding to 10 8 — 10° K) are 
very low, the most relevant excitations in the outer crust are the electrons (specific 
heat c v oc T) and the phonons of the crystal lattice (c v oc T 3 ). In the inner crust, 
the situation is more complicated, since there are in addition the unbound neutrons 
of the gas between the clusters. If these neutrons were a normal Fermi gas, their 
contribution to c v at low T would be linear in T, like that of the electrons but much 
larger due to the higher density of states. However, for most of the relevant densities 
and temperatures one can assume that the neutrons are superfluid (although the 
density dependence of the critical temperature T c of neutron matter is not precisely 
known). In this case, the energy for the creation of a neutron quasiparticle is of the 
order of the pairing gap A ~ 1 MeV, and analogously to the specific heat of a 
superconductor^ the contribution of neutron quasiparticles to the specific heat is 
exponentially suppressed (oc e _A/<T ) at low T. It was therefore argued that whether 
neutron pairing is stronger or weaker can have an observable effect on the cooling 
curve P® Vice versa, observation of neutron-star cooling might help to constrain 
the superfluid critical temperature in neutron star matterP^l P^l 

In contrast to the gapped neutron quasiparticles, long-wavelength collective 
modes of the neutron gas can be thermally excited at low temperature and con¬ 
tribute to the specific heat and the heat conductivity (although the electron and 
lattice-phonon contributions are usually dominant 41 ). 

This article is divided into two parts. In Sec. 0 we start by discussing collective 
modes in homogeneous low-density neutron matter and in particular the appearance 
of the Goldstone mode as a phase oscillation of the superfluid order parameter 
(gap) and its connection with superfluid hydrodynamics. For inhomogeneous crust 
matter, we then summarize results of completely microscopic calculations, which 
are however limited to wavelengths smaller than the distance between neighboring 
clusters, as well as results in the long-wavelength limit. In Sec. 0 we discuss in 
detail a model for collective modes in the pasta phases, in particular in the lasagne 
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phase, based on superfluid hydrodynamics. 

Throughout the article, we use natural units with h = c=kB = l{h = reduced 
Planck constant, c = speed of light, ks = Boltzmann constant). 


2. QRPA, Goldstone Mode, and Superfluid Hydrodynamics 
2.1. QRPA in uniform neutron matter 

For simplicity, let us start our discussion with a uniform neutron gas, although this 
is of course a very incomplete model of the inner crust since it neglects the presence 
of dense clusters. 

Collective modes are conveniently described in the framework of the Random- 
Phase Approximation (RPA) for the response of uniform matter P^J^It has already 
been used in astrophysical contexts, e.g., to calculate the neutrino mean-free pathfSI 
However, the RPA does not include pairing and superfluidity. The extension of the 
RPA to systems with pairing is called Quasiparticle RPA (QRPA). 15 QRPA cal¬ 
culations in uniform neutron and neutron-star matter, using the so-called Landau 
approximation for the particle-hole (ph) interaction, were done in Refs. [T6llT7l Here 
we will discuss results for neutron matter obtained in Ref. [18] with the full Skyrme 
interaction in the ph channel. Also in other fields of physics, the RPA was gen¬ 
eralized to systems with pairing. Let us mention the seminal work by Anderson 19 
and Bogoliubov 2 ® for the case of superconductors, or some more recent work on 
collective modes in ultracold Fermi gases in the unitary limit. 21 22 

The QRPA describes small oscillations around the Hartree-Fock-Bogoliubov 
(HFB) ground statcB and can be derived by linearizing the time-dependent HFB 
equations? 1 ® In the left panel of Fig. [T] we show the QRPA density response func¬ 
tion n(w, q) of neutron matter, as a function of the excitation energy w for different 
values of the momentum q. At first glance, the response resembles the broad ph 
continuum of the usual RPA response. Since the Landau parameter Fq relevant 
for the density channel is negative, there is no zero-sound above the continuum. 23 
However, at low energies, one sees that, because of pairing, the ph continuum start¬ 
ing at uj = 0 is transformed into a two-quasiparticle (2qp) continuum starting at 
a threshold given by approximately 2 AHj This gives rise to the appearance of an 
undamped collective mode below the threshold, indicated by the arrows. 

The nature of this collective mode is a phase oscillation of the gap A = \A\e^. In 
the ground state, the phase is arbitrary but constant (spontaneous 17(1) symmetry 
breaking). Usually it is assumed without loss of generality that A is real, i.e., 
(f> = 0, but a global change of the phase does not cost any energy. However, if the 
phase varies spatially as (f> oc e lq ' r , the excitation energy u> is proportional to q. 
Such a mode related to a spontaneously broken symmetry is called a Goldstone 


a In a uniform gas, the HFB ground state coincides with the Hartree-Fock (HF) + BCS one. 
b Quantitatively, the BCS result for A used in Ref. 1181 is too high since it does not account for 
screening of the pairing interaction! 24 1 2 5 ] 
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Fig. 1. Left panel: QRPA results of Ref. 1181 for the density response II for neutron matter with 
density n = 0.016 fm —3 . The imaginary part of II(a;, q) is shown as function of the excitation energy 
uj for three different values of q. The arrows represent 5-functions corresponding to the undamped 
collective mode. Right panel: Dispersion relation of the collective mode obtained within QRPA 
(solid and dashed lines) and the linear dispersion relation of hydrodynamic sound (dash-dotted 
line). The 2qp thresholc0is indicated by the dots. 


mode,^ 2 '^ in the case of a superfluid it is also known as Bogoliubov-Anderson 
(BA) sound. Note that this BA souncP^Sl exists only in superfluids but not in 
superconductors: In the case of charged particles, the broken symmetry is a local 
one (the electromagnetic gauge symmetry) and in this case there is no Goldstone 
models 

If the phase r) is not constant, the Cooper pairs (of mass 2 m n , where m n 
denotes the neutron mass) move with a collective velocitj^ 


v 





11> • > the Goldstone mode corresponds to a longitudinal density wave. In contrast 
to zero sound, the Goldstone mode does not deform the Fermi sphere during the 
oscillation. Therefore its speed of sound u is given by the hydrodynamic formula 


n n dn n 
m n dn n ’ 


( 2 ) 


where fi n and n n are the neutron chemical potential and density, respectively. 

In the right panel of Fig. [TJ we show the dispersion relation u(q) (solid line) of 
the collective mode obtained in QRPA. As anticipated, it agrees perfectly with the 
hydrodynamic one co = uq for small q (dash-dot line). However, when co approaches 
the 2qp theshold (dotted line), it deviates from the linear dispersion law. At high 
values of q, the collective mode enters the 2qp continuum and gets damped (dashed 
line). 

At low temperatures T < C A, one may neglect the temperature dependence 
of the dispersion relation and the damping of the Goldstone mode. Since mainly 
modes with w ~ T are excited, the deviations from the linear dispersion law ui = uq 
are negligible, too. Therefore, the neutron-gas contribution to the specific heat at 
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Fig. 2. Left panel: HFB density profile calculated by Khan et alP^l in a WS cell with Rws = 
27.6 fm, containing 1750 neutrons and 50 protons ( 1800 Sn). Right panel: quadrupole response 
of neutrons in the same WS cell calculated by Khan et alP*2J within the QRPA (solid line; for 
comparison, the dashed line shows the unperturbed HFB quasiparticle response) as function of 
the excitation energy uj. 


low temperature can be written analytically as 

2tt 2 T 3 


Cy,coll — 


15 u 3 


( 3 ) 


analogous to the specific heat of phonons in a crystal On the one hand, this 
is of course much larger than the contribution of neutron quasiparticles, which 
is suppressed by an exponential factor e~ A ^ T . On the other hand, it is still much 
smaller than the specific heat of normal-fluid (i.e., unpaired) neutrons, which would 
be linear in T. 


2.2. QRPA in a Wigner-Seitz cell 

So far, we have discussed only a homogeneous neutron gas. In order to describe 
clusters in the gas, often the Wigner-Seitz (WS) approximation is employed. This 
approximation consists in replacing the elementary cell of the crystal lattice by a 
sphere of radius Rws, having the same volume as the elementary cell, with the 
cluster in its center. The advantage is that then all quantities, i.e., the densities, 
the mean-field, the Coulomb potential, etc., depend only on the distance r from the 
center of the cell. This spherical symmetry makes it possible to carry out HF 2 or 
HFlpLlLH calculations with the large numbers of particles contained in a WS cell. 
As an example, we display in the left panel of Fig. [2] the HFB density profile in a 
WS cell containing 1750 neutrons and 50 protons ( 1800 Sn) calculated by Khan et 
al.SS The nuclear cluster surrounded by the gas is clearly visible, the gas density 
being practically constant for r > 12 fm (except for a small bump near Rws which 
is an artefact due to the boundary conditions 2 31 ). 

Based on the HFB ground state, Khan et alpZI calculated the collective modes 
within QRPA, too. As an example, the right panel of Fig. [2] shows the quadrupole 
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strength (imaginary part of the quadrupole response function) within the WS cell 
as a function of the excitation energy w. A strongly collective mode, called “super¬ 
giant resonance” in Ref. 1321 appears well below the unperturbed 2qp quadrupole 
excitations (dotted line). 

This super-giant resonance is mainly an excitation of the neutron gas. Its energy 
of ~ 3 MeV as can be seen in Fig. [2] can be roughly understood by considering it 
as the lowest eigenmode of the hydrodynamic (BA) sound in the spherical cell. In 
Ref. (32j this energy was estimated as w = VFq/V3 (vf being the Fermi velocity) 
with q determined from ji(qRws) = 0 ( ji being a spherical Bessel function; in the 
case of the quadrupole mode: l = 2). Note, however, that vf/V 3 is the speed of 
sound of an ideal Fermi gas, which is about 50% higher than the speed of sound u 
one gets from Eq. ©■ 

The argument for the (approximate) validity of hydrodynamics is that the WS 
cell is much larger than the coherence length £ in the neutron gasf 32 Actually, up 
to some factors of order unity, the condition £ <C Rws or £q <C 1 is equivalent to 
ui <C A. Obviously these conditions are not very well satisfied, so that one expects 
to find at a quantitative level some deviations from hydrodynamics, similar to the 
deviation of the QRPA dispersion relation ui(q) in uniform matter from the linear 
one, u) = uq (cf. Fig. []]). 

Let us mention that the correspondence between QRPA and superfluid hydro¬ 
dynamics in non-uniform systems in the case of strong enough pairing was demon¬ 
strated in Ref. [331 in the context of trapped ultracold Fermi gases, too. However, 
because of the typically very large numbers of atoms and strong pairing in these 
systems, the situation is usually much more favorable than it is in the neutron star 
crust, and superfluid hydrodynamics can give precise predictions for the frequencies 
of collective modesP^ 


2.3. Low-energy theory for large wavelengths 

While the WS approximation is well suited for the description of static properties, it 
is obviously not capable of describing collective modes whose wavelengths exceed the 
size of the WS cell. In reality, however, the most relevant modes at low temperature 
are acoustic modes, i.e., modes whose energy oj is proportional to q for q -C 1 /Rws- 
As mentioned in Sec. 12.11 the Goldstone mode in the uniform neutron gas is a 
consequence of the broken global U{ 1) symmetry and corresponds to an oscillation 
of the phase (j> of the neutron gap. This argument remains valid in the presence of 
clusters. This lead Cirigilano et alP^to develop an effective theory for the combined 
system of the superfluid neutron gas and the crystal lattice of clusters. The degrees 
of freedom, i.e., the fields appearing in the effective Lagrangian, are the phase <j>( r, t ) 
of the neutron gap and the displacement £(r .t) of the clusters. Similar low-energy 
theories were subsequently used in Refs. |39]!40] Since <j> and £ are coarse-grained 
over regions larger than the periodicity L of the crystal lattice, these low-energy 
theories are only valid for q <C 1/L. 
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Fig. 3. Velocities of the longitudinal and superfluid phonons (m and u^, left panel) and of the 
transverse phonons ( ut , right panel) as a function of the baryon density ng, from Charnel et alP^ 
Dotted lines: entrainment neglected; dashed lines: entrainment included; solid lines: entrainment 
and mixing included. 
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with entrainment: u t 


•r v 



Of course, clusters and neutron superfluid do not move independently of each 
other. This is called the entrainment effect.^ It should be noted that in contrast to 
an ordinary dragging of the gas by the clusters (and vice versa), the entrainment 
is a non-dissipative force. It modifies the velocities ui and u t of the longitudinal 
and transverse lattice phonons, respectively, as well as the velocity u</, of the Gold- 
stone mode (also called superfluid phonon) Furthermore, the longitudinal lattice 
phonons and the Goldstone mode get mixed 38 ® and one finds two new eigenmodes 
with two new velocities u ±. The framework of low-energy effective theory was not 
only used to describe the long-wavelength phonons, but also phonon-phonon and 
phonon-electron interactions. These are particularly relevant for the calculation of 
the phonon (lattice and superfluid) contribution to transport properties such as the 
heat conductivity.®] 

The coefficients of the effective theory, in particular the entrainment and mixing 
coefficients, were calculated by Charnel 37 ^ using a band-structure calculation for the 
neutrons in the periodic mean field of the cluster:!^ analogous to those used in solid- 
state physics for the electrons in the periodic Coulomb field of the ions in a metal. 
As an example, we display in Fig.[3]the results of Ref. [39jfor the density dependence 
of the different velocities without (dotted lines) and with (dashed lines) entrainment 
and with both entrainment and mixing (solid lines). One observes that the coupling 
between the lattice and the superfluid neutron gas has a drastic effect on the lattice 
phonons. The entrainment results in a reduction of the phonon velocities, since some 
neutrons of the gas move together with the clusters, increasing the cluster effective 
mass (in the language of Refs. 137(1391 the density of “conduction neutrons” moving 
independently of the clusters is smaller than the density of “unbound” neutrons). 

However, as pointed out by Kobyakov and Pethick®! zero-point oscillations of 
the clusters may reduce the band structure effects. Furthermore, we note that, using 
a completely different approach (similar to the one we will use in the next section), 
Magierski and BulgacP^® arrived at the opposite result, namely that the effective 












June 2, 2015 0:46 WSPC/INSTRUCTION FILE paper 


8 M. Urban and M. Oertel 



Fig. 4. Periodic slab structure as simplified model for the lasagne phase. 


mass of the clusters is reduced when they are immersed in the neutron gas. We 
conclude that the microscopic modeling of the coefficients of the effective theory, 
in particular those depending on the entrainment, is not yet completely settled. 


3. Hydrodynamic Model for Collective Modes in the Pasta Phases 

After this survey of works pointing out the connection between the Goldstone mode 
of the superfluid neutrons and superfluid hydrodynamics, let us now describe in 
some detail our hydroynamic model for the collective modes, extending our preced¬ 
ing work. 4, 3ln particular, we include the Coulomb interaction, neglected in Ref. [451 
in order to describe simultaneously superfluid modes and lattice phonons. Our aim 
is to cover also wavelengths that are large compared to the coherence length (cf. 
discussion in Sec. 12.211 but not necessarily large compared to the periodicity L of the 
crystal lattice. In this sense the approach bridges between the two extreme cases 
discussed above. In principle the formalism can be applied to all phases of the inner 
crust, but in the practical calculations we restrict ourselves to the simplest case, 
which is the phase of plates (“lasagne”), see Sec. 13.41 


3.1. Nuclear pasta as phase coexistence in equilibrium 

For the ground-state configuration, we use a very simple model. We describe matter 
in the inner crust as a mixed phase consisting of a neutron gas (phase 1) and a 
neutron-proton liquid (phase 2). We assume that in equilibrium the densities in each 
phase are constant and the two phases are separated by a sharp interface. We neglect 
the smooth variation of the density, found in microscopic calculations (cf. left panel 
of Fig. |2]), at the transition from one phase to the other. Since the phases coexist, 
the pressure and the chemical potentials must be equal in both phases. In addition, 
an electron gas globally compensates the charge of the protons. We approximate it 
by an ideal gas of massless fermions with uniform density. Since neutrons, protons 
and electrons are in /3-equilibrium, their chemical potentials satisfy fi n i = p n 2 = 
hp 2 + he ■ Obviously the model as stated above is very crude. In order to account at 
least in some approximate way for the realistic (smooth) interface, we assume that 
its effect can be subsumed in a single parameter, the surface tension a (energy per 
interface area). 
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3.2. Hydrodynamic equations with Coulomb potential 

Since neutrons and protons are paired, we assume that their motion can be de¬ 
scribed by superfluid hydrodynamic equations, except at the interface between the 
two phases, where appropriate boundary conditions are needed. This is similar to 
the approach used in Refs. EE51PH1 to calculate the flow of the neutron gas around 
and through a moving cluster and the corresponding effective mass of the cluster. 

In superfluid hydrodynamics, the momentum per particle of each fluid, p a (a = 
n,p indicating neutrons or protons), is related to the phase </> a of the corresponding 
superfluid order parameter by 

Pa(r,t) = . (4) 

In the neutron gas, the momentum p„ is proportional to the velocity v n and 
we may write 

P n = M n V n . (5) 

In the non-relativistic limit, i.e., if the total energy density is dominated by the 
mass density, the proportionality constant M n is of course given by the neutron 
mass m n , and one recovers the well-known expression m for the superfluid velocity. 
Relativistically, however, one findJ^that M n is equal to the enthalpy per particle, 
which, at zero temperature, equals p, n . 

In the liquid phase containing neutrons and protons, the situation is more 
complicated. If both fluids have different velocities v n ^ v p , the interaction be¬ 
tween neutrons and protons can give rise to entrainment, i.e., to a misalignment 
between momenta and velocities,^ even in a uniform system and already at a mi¬ 
croscopic level (i.e., independently of the macroscopic entrainment effect mentioned 
in Sec. 12.31 which is a consequence of the crystalline structure). This microscopic 
entrainment effect can easily be understood in the framework of Landau Fermi- 
liquid theoryP*3IB>® q^e relationship between velocities and momenta can then be 
written as 

Pa = ^ M abVb , (6) 

b=n,p 

where the matrix M is related to the so-called Andreev-Bashkin entrainment ma¬ 
trix, Y in the notation of Ref. m via M a b = Yhb^ab n b- I n the non-relativistic 
case, Galilean invariance implies that Y2b M a b ~ w °’ an( I the absence of entrain¬ 
ment, the matrix M reduces to M a b = m a 5 a b■ In the relativistic case, M satisfies 
Y2b M ab = ha (see I Appendix A[ ). 

After these preliminary remarks, let us turn to the hydrodynamic description 
of small oscillations. The first hydrodynamic equation is the continuity equation, 
describing the conservation of particle number of each species a. Since we are only 
interested in small oscillations, we consider the fluid velocities v Q as small. In ad¬ 
dition, we decompose the density into its constant equilibrium value n a , e q and a 
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small deviation Sn a (r,t) = n a (r,f) — n a , e q- Keeping only terms linear in the small 
deviations from equilibrium, we can write the continuity equation for species a as 

d t Sn a {r, t) = -n a , eq V • v a (r, t). (7) 

For the sake of better readability, we will from now on drop the index eq since it is 
clear that, after linearization, any quantity that is multiplied by a small quantity 
(such as v a or Sn a ) has to be replaced by its equilibrium value. 

The second equation is the Euler equation, describing the conservation of mo¬ 
mentum. Keeping again only terms linear in the deviations from equilibrium, we 
obtain in the case of the neutron gas 


dtp n (r,t) = ~V<5n„ , 
on n 


( 8 ) 


where the derivative has to be taken at the equilibrium density n n . In the liquid 
phase, neutrons and protons are coupled by strong interactions, and in addition the 
protons feel an acceleration due to the variation of the Coulomb potential, SV. The 
corresponding two (a = n,p) Euler equations are 

dfJ-a 


dn b 


VSn b (r, t) - SapVSV (r, t). 


( 9 ) 


dtp a {r,t) = - Y 

b=n,p 

Note that the relations between v a and p a , Eqs. © and ©, are linearized, too, 
i.e., it is sufficient to calculate M at the equilibrium densities. 

The variation of the Coulomb potential depends itself on the variation of the 
proton and electron densities. We assume that the electron density follows instanta¬ 
neously the motion of the protons, leading to a screening of the Coulomb interaction. 
The corresponding modified Poisson equation for the Coulomb potential reads 


(V 2 - -i^<5K(r, t) = -4ite 2 5np 


( 10 ) 


(ID 


where A is the Debye screening length given by 

1 , 0 dn e 

w = ne dUe ■ 

Let us now consider harmonic oscillations, i.e., all deviations from equilibrium 
(<t> a , p a , v a , Sn a , SV) oscillate like e~ lult . Therefore we can replace all time deriva¬ 
tives d t by a factor — iw. Using Eqs. © and 0, we can express v Q and Sn a in 
terms of <p a and obtain the following equations for (f> a and SV : 

(i) In the gas phase without protons: 


V 2 ^, =- 


2 'rn i 


1 


V 2 SV=-^SV, 


( 12 ) 

(13) 


where 


u 2 = 


n n dpr, 
M n c)tIji 


(14) 
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denotes the square of the sound velocity, 
(ii) In the liquid phase with protons (a = n,p ): 


= - W 2 ]T Ujct> b -2iuU- p 2 5V, 


b—n,p 


b=n,p 0 ^ 


where we have defined the abbreviations 

1 dn c 
n c 5 /if, ’ 


U ab = £ M * 


c=n,p 


— = 4ne 2 ^ 
A 2 dfj, a ■ 


(15) 

(16) 


(17) 


Note that the eigenvalues of the matrix U are the sound velocities of the two 
eigenmodes in the uniform neutron-proton liquid without Coulomb interac¬ 
tion!^ 


3.3. Boundary conditions 

As mentioned before, the equations of the preceding subsection are valid inside each 
phase but not at the interface between the gas and the liquid phase. At the inter¬ 
face, they have to be supplemented by suitable boundary conditions. In Ref. (45j we 
used very simple boundary conditions: the pressure and the normal velocities had 
to be equal on both sides of the interface. These conditions correspond to an imper¬ 
meable interface. In the present work, we will improve the boundary conditions and 
allow for a neutron flux across the interface, as in Refs. [4310 Similar boundary 
conditions were also given in Ref. [48] in a completely different context, namely for 
an unpolarized Fermi gas surrounded by a polarized one in an atom trap. 

Let us consider the case of a surface parallel to the xy plane, separating the 
gas phase (1) from the liquid phase (2). Since there are no protons in phase 1, the 
velocity of the surface is obviously equal to the normal component of the proton 
velocity in phase 2, v zp 2 , and the displacement of the surface from its equilibrium 
position is given by £ z = v zp /(—iui). If the normal component of the neutron velocity 
is different from the velocity of the surface, this means that neutrons cross the 
surface and pass from one phase to the other. The requirement that the neutron 
current leaving phase 1, n n i(v zn \ — v zp 2), must be equal to the neutron current 
entering phase 2, n n 2(v zp 2 — v zn 2 ), gives our first boundary condition: 

‘fl'nl'Vznl = ‘fl / n2'Vzn2 j (1$) 

where 


— n n 2 n n i . 


This equation can be rewritten in terms of the phases <j) a as 


Ttnl 


d z (t>, 


n 1 


(' n n2 M n * - &n n M p a)d z (t) a 2 . 

a=n,p 


(19) 


( 20 ) 
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In our notation, the quantity M n belongs to the neutron gas and is calculated at 
density n n i, while M a f, belongs to the liquid phase and is calculated at densities 
n n 2 and n p2 - 

Furthermore, the fact that neutrons can cross the surface implies that, as in 
equilibrium, the neutron chemical potentials on both sides of the surface must be 
equal, i.e., 

= SfJ, n 2 . ( 21 ) 

If we express <5/r„ in terms of derivatives d^ a /drib and the density variations Sni, 
and use in each phase the equations of Sec. 13.21 Eq. (12U) reduces to the very simple 
condition 


4>nl — <^n2 ■ (22) 

The condition Pi = P 2 of equal pressure^] used in Ref. [¥5] becomes more com¬ 
plicated once the surface tension is included. The pressure difference is then given 
by the Young-Laplace formulaP® 

P2 ~ Pl=a (i^ + i^)’ (23) 

where Pi and P 2 are the principal curvature radii of the interface and the sign 
is such that the pressure is higher in the phase having a convex surface. Since we 
consider here a surface that is flat in equilibrium, its curvature arises only from the 
displacement £, and can be expressed in terms of derivatives of 4> a2 - The pressure 
difference on the left-hand side of Eq. (l23l) is related to the density oscillations, 
for instance we can write the deviation of the pressure in the neutron gas from its 
equilibrium value as 6P\ = n n iSfi n i. Using again the equations of Sec. 13.21 one can 
rewrite Eq. as 

w( ^2 n “2</>a2 - «n!0ni) + 2 in p2 SV2 = ±— ^2 M~^(d% + Oy)d z (j ) a 2 , (24) 

a=n,p a=n,p 

where the upper (lower) sign is valid in the case that phase 2 is situated above 
(below) phase 1. 

Let us now turn to the boundary conditions related to the Coulomb potential 
V. First of all, V itself has to be continuous at the surface, i.e., 

SVl = SV 2 . (25) 

This ensures that the component of the electric field tangential to the surface is 
continuous, too. 47 

To linear order in the deviations, the charge of the protons in the region between 
the unperturbed and the perturbed surface can be considered as a surface charge 


c To linear order in the velocities, the distinction between pressure and generalized pressurJ^H^l 
is irrelevant. 
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density en P 2 ^. This gives rise to a discontinuity of the electric field normal to the 
surface,^ 

SE Z i - 5E Z2 = 47ren P 2^2 ■ (26) 

Expressing again all quantities in terms of the potentials 5V and </> Q , we obtain our 
last boundary condition 

-2iojd z (5V 2 -SV 1 )= n 2 a d z ct>a 2, (27) 

a=n,p 

with the abbreviation 

fl 2 a = 4tt e 2 n p2 M~ c ) . (28) 


3.4. Collective modes in a periodic slab structure (lasagne phase) 

We will now consider the simple case of a periodic structure of slabs. We suppose 
that for 0 < z < L\ we are in the gas phase (1), for Li < z < L = L\ + L 2 we are 
in the liquid phase (2), and for L < z < L + L\ we are again in the gas phase, and 
so on, as shown in Fig. 0] 

Since this system is translationally invariant in x and y directions, it is clear 
that the x and y dependence of the eigenmodes (i.e., of the potentials cj> a and SV) 
is of the form 

</> a (r) = ( f )a {z)e i ^ x+q y y) , SV(r) = 5V(z)e i{q * x+q * y) . (29) 

From the periodicity of the system in z direction it follows that the eigenmodes 
satisfy the Bloch conditional 

<t>a(z + L) = (j>a{z)e iq * L , SV{z + L) =SV(z)e iq * L . (30) 

For given q, these conditions can only be satisfied for some discrete values of u. This 
is the excitation spectrum we are looking for. Because of the Bloch conditions (l30|) . 
it is sufficient to solve the coupled differential equations of Sec. 13.21 in the regions 
(1) from 0 to L 1 and (2) from Li to L, together with the boundary conditions of 
Sec. 13.31 at, z = L\ and L. 


3.5. Numerical results 

Let us now investigate the resulting excitation spectrum for a specific example. 
The values for the equilibrium quantities will be taken from the work by Avancini 
et al.p who have studied the structure of pasta phases in a relativistic mean 
field model. Our geometry corresponds to the lasagne phase, which has been found 
in Ref. [42j in the case of zero temperature and /3-equilibrium for baryon number 
densities 0.077 fm -3 < r 13 0.084 fm -3 . For our example we have chosen an 

intermediate density, ns = 0.08 fm” 3 , as in Ref. [45) The corresponding properties 
of the two phases (1) and (2) are summarized in Table [U 
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Table 1. Properties of the lasagne phase within the model by Avancini et 
al.^1 studied in our example. The average densities of the total system are 
given by n n = ( L\/L)n + {L 2 /L)nffi etc. Baryon density and proton 
fraction are defined as ng = n n + n p and Y p = n p /n n , respectively. 





slab (1) 

slab (2) 

total 


L 

(fm) 

9.40 

7.38 

16.78 


Tin 

(fm- 3 ) 

0.0701 

0.0885 

0.0782 


n p 

(fm” 3 ) 

0 

0.0041 

0.0018 

n B 

—— n-fi + n p 

(fm” 3 ) 

0.0701 

0.0926 

0.0800 

Y P 

= n p /n B 


0 

0.0447 

0.0227 


The microscopic results for the equilibrium configuration will be used to estimate 
the value of the surface tension, too. While the surface tension favors big structures, 
the Coulomb interaction favors small ones. In the case of the Lasagne geometry, 
the Coulomb and surface energies per volume are given by 

ec + es = -^e 2 n 2 e L\ + — , (31) 

where n e is the electron number density, related to the proton number density in 
phase 2, n p 2 , by the requirement of charge neutrality ( n e L = rip^L'i). and L = 
L\ + L 2 . The size L of the structure can be obtained by minimizing ec + es with 
respect to L, keeping the ratio L 1 /L 2 fixed. Inversely, the value of the surface 
tension can be estimated from the values of L± and L 2 given by the microscopic 

calculations P as 

<7 = L 2 rig L L\ . (32) 

6 

In principle the Coulomb interaction will result in slightly non-uniform density 
distributions of the protons inside phase 2 and of the electron gas. However, from 
the Coulomb potential calculated with the uniform density distributions one can 
estimate that this effect is weak and can be neglected in a first approximation. 

The same remarks of caution as in Ref. [45] concerning the dimensions of the 
structure, the coherence length and the validity of the superfluid hydrodynamics 
approach (cf. discussion in Sec. 12.211 are of course valid here and the interested 
reader is refered to that paper. 

Let us now discuss the solutions for the energies oj and compare the present 
model with that of Ref. 051 In Fig. [5] the energies are shown as functions of 
q for three different angles 9 between q and the 2 axis (i.e., q z = qcosd and 
911 = + q p = qsm9). The left panel shows the dispersion relation for waves 

propagating in z-direction, i.e. perpendicular to the interfaces between the different 
slabs. We observe an acoustic branch with an approximately linear dispersion law 


ui = Uiq 


( 33 ) 
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8 = 0 ° 


8 = 45 ° 


6 = 90 ° 




q (fm' 1 ) q (fm 4 ) q (fm 4 ) 


Fig. 5. Dispersion relations of the modes propagating along the 2 -axis (9 = 0, left), at an angle of 
45°, and in the plane parallel to the slabs (9 = 90°, right). The plain red lines correspond to the 
present model and the dashed blue lines to the model of Ref. [451 neglecting the Coulomb interaction 
and assuming an impermeable interface between the different phases. The dotted magenta lines 
indicate the results of the present model upon neglecting the microscopic entrainment effect and 
the dash-dotted turquoise line the result neglecting the surface tension. 


at low energies, and several optical branches with a finite energy for q = 0, analo¬ 
gously to phonon branches in a crystal. The main differences between the present 
model (solid red lines) and the previous onc^H (dashed blue lines) are as follows: 

(i) Here we use more realistic boundary conditions corresponding to a permeable 
interface between the two phases, whereas in Ref. |45] neutrons were not allowed 
to pass from one phase into the other. 

(ii) We include the Coulomb interaction which was neglected in Ref. [45, 

Actually, the change of the boundary conditions has only little influence on the 
energy spectrum, whereas the Coulomb interaction gives rise to an additional mode. 
Compared with Ref. 051 we have introduced the surface tension, too. The effect of 
the surface tension is obviously vanishing for modes propagating perpendicular to 
the interface and maximal for modes propagating along the layers. The effect on the 
excitation energies is small, see the right panel of Fig.0 where we display the result 
for a = 0 in comparison with the complete calculation. On the other panels the 
difference is too small to be seen and therefore not shown. The effect of microscopic 
entrainment is weak, too, as can be seen from the left panel. 
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Fig. 6. Specific heat as a function of temperature for ns = 0.08 fm 3 . 


Note that within the Wigner-Seitz approximation as discussed in Sec. 12.21 we 
would only obtain a discrete spectrum corresponding to our spectrum in the case 
<7 = 0. The reason is that in this approximation the coupling between cells is 
neglected, and thus each cell has the same excitation spectrum. The degeneracy of 
the modes in each cell is lifted by the coupling between cells, which gives rise to a 
momentum dependent spectrum as obtained in our approach. 

Let us now look at the cases 9 = 45° and 90° shown in the central and right 
panels of Fig. [5l As in Ref. |45j we find now two acoustic modes. The energy of the 
second acoustic mode can be approximately written as 

w = u 2 q\\ = u 2 gsin0, (34) 

which means that this mode propagates in a slab with almost no coupling to neigh¬ 
boring slabs. 

The two acoustic modes dominate largely the baryonic contributions to the 
specific heat and exceed in particular the contributions from superfluid neutron 
quasiparticles. The temperature dependence of the specific heat corresponding to 
the modes shown in Fig. [5] is displayed in Fig. [6] We assume here that the temper¬ 
ature dependence of the mode spectrum can be neglected, which should be true for 
T <C A„, A p . The electronic contribution c v>e = /z 2 T/3 and the contribution from 
weakly paired neutron quasiparticleJ® are shown for comparison, too. 

At very low temperatures, the specific heat of the collective excitations is ap¬ 
proximately given bj® 

c V: coii ~ aT 2 + bT 3 , (35) 

with a ss 3C(3)/7 tu|L and b ~ 27r 2 /15uf. The T 3 term is the usual acoustic phonon, 
whereas the “two-dimensional” mode discussed above gives rise to a contribution 
increasing as T 2 PS We note, however, that within our new model the deviations 
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from this approximate formula are stronger than in Ref.[45]for both acoustic modes: 
the level repulsion due to avoided crossings let them deviate from the purely linear 
behavior for much smaller momentum values. As a consequence, the approximation 
(135[) is valid only for much smaller temperatures, where the electronic contribution 
to the specific heat is the most important one. As can be seen from Fig. [6l the 
additional mode and the resulting flattening of the acoustic modes increase the 
contribution to the specific heat. 


4. Summary and Conclusions 

In Sec. El we reviewed a couple of existing approaches to describe collective modes 
in the superfluid inner crust of neutron stars. Considering only a uniform neutron 
gas within the QRPA, one recovers the Goldstone mode, which is a longitudinal 
density wave (BA sound) with linear dispersion relation at small q , as the lowest 
lying collective excitation. For a realistic description of the neutron star inner crust, 
however, its crystalline structure should be considered, i.e. it should be described as 
a lattice of dense clusters surrounded by a neutron gas. There are several approaches 
to treat collective excitations of this non-uniform matter involving different kinds 
of approximations. 

In the standard hydrodynamic approach, only coarse-grained quantities are con¬ 
sidered, i.e., quantities averaged over many cells. This leads to an effective low- 
energy theory, in which the superfluid collective modes with linear dispersion re¬ 
lation are coupled to the dynamics of the crystal lattice. However, this effective 
theory is only valid at very large wavelengths (much larger than the periodicity of 
the lattice, L), and the determination of its coefficients, especially those involving 
the so-called entrainment, is difficult. 

On the other hand, a QRPA calculation within the WS approximation resolves 
the microscopic structure of each cell. However, it gives the lowest lying collective 
mode at a finite energy. This is an artefact of treating an isolated cell with a finite 
radius Rws- Nevertheless the result is very important since it shows that in a non- 
uniform system, superfluid hydrodynamics may be used at a microscopic level (i.e., 
on length scales smaller than Rws) as long as the coherence length is sufficiently 
short, i.e., pairing is sufficiently strong. 

In Sec. [3] we have presented the approach of superfluid hydrodynamics at a 
microscopic level in order to develop a model for the collective modes of large and 
intermediate wavelengths (i.e., larger than or comparable with L ). In this model, 
hydrodynamic equations in the clusters and in the neutron gas were matched at 
the phase boundaries with boundary conditions allowing neutrons to pass from one 
phase into the other. For simplicity, the model has so far only been solved in the 
phase of plates. The most striking result is the appearance of an approximately 
two-dimensional collective mode, which gives a contribution to the specific heat 
proportional to T 2 (in contrast to the usual T 3 behavior) at low temperatures. 

There remain many open questions. In particular, in order to better understand 
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the two-dimensional mode, we are developing an effective theory corresponding 
to our model in the limit of small q. This could help to clarify some questions 
concerning entrainment, too see the discussion in Sec. 12.31 In addition, we 
plan to apply our hydrodynamic model to other geometries such as the crystalline 
phase or the phase of rods. 

The latter point is essential in order to make contact with astrophysical ob¬ 
servations, since it would allow one to determine heat transport properties in the 
entire inner crust, necessary for modeling neutron star thermal evolution. Proper¬ 
ties of the crust influence the cooling process mainly during the first 50-100 years, 
when the crust stays hotter that the core which cools down very efficiently via neu¬ 
trino emission (see, e.g., Ref. [SJ. Heat transport in the crust is the key ingredient 
to explain the afterburst relaxation in X-ray transients, too.® ^ As microscopic in¬ 
gredients for the models of the thermal relaxation of the crust, in addition to the 
specific heat, thermal conductivity would be needed, and to less extent neutrino 
emissivities. Although the contribution of collective modes to heat transport prop¬ 
erties exceeds largely that of superfluid quasiparticles, it is to be expected that in 
almost the whole crust electronic contributions remain dominant, except in pres¬ 
ence of a neutron star magnetic field of the order 10 13 G or higherpS which is the 
case for many observed neutron stars. Therefore it would be interesting to extend 
the present model to a magnetised environment, too. 


Acknowledgements 

We thank Elias Khan for providing us with the numerical results of Ref. 1521 


Appendix A. Microscopic Input for the Model of Section [3] 


As microscopic input, we need the equation of state, i.e., the relation between the 
densities n a and the chemical potentials p a . As in Ref. 1451 in our concrete numerical 
examples, we use the results of the work by Avancini et al. 42 for the equilibrium 
configurations. They evaluate the structure of the pasta phases for charge neutral 
matter in /3 equilibrium using a density dependent relativistic mean-field model, 
the DDH<5 model (originally called DDHpJ), for the nuclear interactionP^^Hl i n 
order to be consistent, we shall calculate the chemical potentials p a with the same 
interaction. The entrainment coefficients within this model, evaluated following the 
approach in Refs. [SHIM} are given by 


t \t _ Lr i f^R . I u , t a tb i p 
Mab — RT'a ®ab + Rb I-1-T H- Z -y 

V ub mf, 4 mi 


t a t b r I 


(A.l) 


where ub = n n + n p denotes the baryon density, t n = — 1 and t p = 1. The Landau 
effective masses are denoted = \fp 2 Fa + not to be confused with the Dirac 
effective masses rn p . The Fermi momenta are represented by p Fa and Tj = gihi{ns) 
are the meson-nucleon coupling constants of the DDHd model, see Ref. [42] The term 
Y,r, not present in the expressions given in Refs. 151(1521 arises from the density 
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dependence of the coupling constants and appears in the Landau effective masses, 
too, 

r 2 r 2 

nia = Fa - ^ n B ~ t a —^( n p ~ n n ) - Y, R . (A.2) 

m u 4TO p 

and an explicit expression can be found in Ref. (421 
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